Noise-induced clustering in Hamiltonian systems 



o 
o 



D.V. Makarov, M.Yu. Uleysky, M.V Budyansky, and S.V. Prants 
Laboratory of Nonlinear Dynamical Systems, 
V.I. Il'ichev Pacific Oceanological Institute of the Russian Academy of Sciences, 
690041 Vladivostok, Russia 
(Dated: February 8, 2008) 

The motion of oscillatory-like nonlinear Hamiltonian systems, driven by a weak noise, is consid- 
ered. A general method to find regions of stability in the phase space of a randomly-driven system, 
based on a specific Poincare map, is proposed and justified. Physical manifestations of these regions 
of stability, the so-called coherent clusters, are demonstrated with two models in ocean physics. We 
find bunches of sound rays propagating coherently in an underwater waveguide through a randomly 
fiuctuating ocean at long distances. We find clusters of passive particles to be advected coherently 
by a random two-dimensional fiow modelling mixing around a topographic eddy in the ocean. 
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I. INTRODUCTION 

The interplay between dynamical chaos and noise is 
a topic of interest both in nonlinear dynamics and sta- 
tistical physics. Real-world systems operate against a 
noisy background. Both chaos and noise mean random 
or diffusive-like behavior in the nonlinear system's dy- 
namics that can be quantified by the maximal Lyapunov 
exponents 



A = lim A(t), X{t) 



lim lln#l. (1) 
A(o)^o t A(0) ' ^ ' 



where A{t) is a distance (in the Euclidean sense) in a 
direction at the moment t between two initially closed 
trajectories. 

However, the methods and approaches to study deter- 
ministic and random system are different. Typical deter- 
ministic nonlinear systems (the so-called nonhyperbolic 
system) are known to have mixed phase space with coex- 
isting regions of stable and unstable motion. There are 
islands of regular motion merged into a chaotic sea, is- 
lands around islands, stochastic layers confined between 
invariant tori, cantory nearby the boundaries of islands, 
and so on (for a recent review of chaos in Hamiltonian 
systems see [J). The motion is predictable in some re- 
gions but permits a statistical description only in those 
ones where extremal sensitivity to initial conditions takes 
place. The notion of the mixing time, the reciprocal max- 
imal Lyapunov exponent, can be introduced to quantity 
the so-called predictability horizon. 

To describe randomly-driven systems one uses some 
ergodic postulates based on the assumption of fully de- 
veloped chaos. In Hamiltonian nonintegrable systems 
intermittent-like dynamics with chaotic oscillations in- 
terrupted by regular ones seems to be more realistic. 
It means, in particular, that regular motion docs not 
cease suddenly, and completely stochastic motion does 
not arise suddenly just beyond the predictability hori- 
zon. Remnants of stability should persist for more long 
times. It is especially true for oscillatory-like systems 



due to resonant interaction of unperturbed motion with 
random perturbation 2\. 

In this paper we treat a model of a randomly-driven 
nonlinear oscillator from purely deterministic point of 
view, and show that there are regions of stability sur- 
viving for a comparatively long time under a weak noisy 
excitation with arbitrary spectrum. This time can be es- 
timated with the help of a specific map which plays the 
role of the Poincare map in periodically-driven systems. 
The map is intended to find numerically clusters of tra- 
jectories with close initial conditions which are stable by 
Lyapunov with a given realization of the random pertur- 
bation (Sec. II). We check the effectiveness of the proce- 
dure proposed with two problems of ocean physics where, 
as in other geophysical problems, there is, as a rule, no 
respective statistical ensemble of averaging, and we deal 
with single realizations in the field experiments. Noise- 
induced coherent clusters are demonstated with the ray 
model of sound propagation in an ocean waveguide with a 
randomly fluctuating sound-speed profile induced by in- 
ternal waves in the deep ocean (Sec. III). Another exam- 
ple is clustering of passive particles advected by a random 
two-dimensional velocity field in an ideal fluid (Sec. IV). 



II. EFFECTIVE POINCARE MAP 

In this section we introduce an effective Poincare map 
that, like to the usual Poincare map, enables to find 
numerically regions of stability in the phase space of a 
Hamiltonian system surviving under a weak random per- 
turbation due to nonlinear resonances between the unper- 
turbed motion and the perturbation. 

Let us consider a one-dimensional nonlinear oscillator 
with the Hamoltonian 

H = Ho + sHi (t) = y + U{q) + eV{q)m , (2) 

where q and p are position and momentum, U(q) is an 
unperturbed potential, V{c[) is a smooth function, S^it) 
is a weak noise with e <C 1, and normalized first and 



2 



second moments, (^) = and = 1/2. Hereafter 

we will consider a single realization of noise, ^(t), and 
therefore, the equations of motion 

dq dp dU dV 



dt ^' dt dq dq 



(3) 



can be treated as deterministic ones. Introducing the 
canonical transformation from the variables to the 

action-angle variables (see P, Q ) 



I=-^j)pdq, ■& ^ j p{q, H) dq, 

we rewrite Eqs.© in the form 

dl _ dH(I, d) _ dV{I, d) 



(4) 



dt dd 
dd _ dH{I, I?) 



dt dl 



(5) 



where uj{I) = dHo/dl is an action-dependent character- 
istic frequency of oscillations. 

The notion of stability means a weak sensitivity of a so- 
lution to small changes in initial conditions that is quan- 
tified by the Lyapunov exponent Q. As we know from 
KAM-theory Q , most of the invariant tori of integrable 
Hamiltonian systems, where the motion is stable, are pre- 
served under a small perturbation. In spite of each single 
realization of a random perturbation ^(t) can be treated 
as a deterministic function than a stochastic one, infi- 
nite number of frequencies in the spectrum of ^(t) leads 
to densely distributed resonances in the phase space and 
laking of invariant tori. In the limit t — > oo, no regions 
in the phase space remain forbidden under a noisy per- 
turbation. Even a weak noise can accelerate phase-space 
transport by forcing trajectories to traverse KAM tori 

IE IS ' So deterministic description ceases to have a 
sense, and one forces to use statistical description of the 
motion. 

However, in practics, we deal with a finite time inter- 
val [0 : Tq]. Moreover, in geophysical field experiments 
experimentalists work mainly with single realizations of 
nonstationary processes of interest. By these reasons we 
can recover "deterministic" methods to explore the sets 
of stable trajectories satisfying to the condition of the 
finite-time invariance: if any set in the phase space at 
t = transforms to itself at t = Tq without mixing, then 
it corresponds to an ensemble of trajectories which are 
stable by Lyapunov within the interval [0 : Tq]. In or- 
der to find such stable sets for an arbitrary spectrum of 
perturbation, we propose the following map 

h+i = I{t - To, = §{t = To, 

(6) 

where I{t; and d{t] li^di) are the solutions of 

Eqs.© with initial conditions /(O) — li, i?(0) = In 
fact, this map is equivalent to a Poincare map for a sys- 
tem with the Hamiltonian 



(7) 



where ^{t) is a periodic function consisting of identical 
pieces of £,{t) of the same duration Tq 



at + nTo)^at), 0<t<To, 



0,1, 



(8) 



In this way we replace our original randomly-driven sys- 
tem by a periodically-driven one. It should be empha- 
sized that the validity of this replacement in restricted 
by the time interval [0 : To]. By analogy with the usual 
Poincare map, the key property of the map © is the fol- 
lowing: each point of a continuous closed trajectory of the 
map 0) corresponds to a starting point of the solution of 
Eqs. |3|j which remains stable by Lyapunov till the time 
Tq. The inverse statement is not, in general, true. It will 
be shown below that the map ^ provides sufficient but 
not necessary criterion of stability. Topological proper- 
ties of trajectories of the map © can be treated in the 
framework of the theory of nonlinear resonance The 
functions V{L, d) and ^{t) can be decomposed in Fourier 
series 



1— — 00 

oo 



(9) 



m— — oo 



where Q ~ 27^ /Tq. The Fourier amplitudes of the analytic 
function V decay as Vi{I) ^ exp(— cri), whereas the ones 



for the random function ^(t) decay as ^„ 



(0< 



/3 < 1 and /? = for a white noise) since ^(0) ^ S,{Tq). 
If the interval Tq is large enough, the parameter (3 is 
determined mainly by the correlation time of S,{t). 

Substituting the expansions © into the equations 
of motion Q and omitting rapidly oscillating terms 
exp[i(M -f mM)], we get the following equations: 



JT - CO 

dl le 

— = /Vi4„i cxpi$ 4- c. c. 



dd e ^ dV, 



(10) 



i,m=0 



where $ ~ W — mil.t+(l)i — ipm- The stationary phase con- 
dition d^/dt = implies the resonances of the map © 



2T{I) - ;To 



(11) 



where T(/) — 2tt lu){l) is a period of unperturbed oscilla- 
tions with a given value of the action /. Resonant values 
of the action /j-cs can be found from the condition l|ll|l . 
The relation I : to defines the order of the respective reso- 
nance. It should be noted that an infinite number of res- 
onances jl : jm (j = 1, 2, ...) corresponds simultaneously 
to each value of the resonant action. However, if /les is 
far enough from the separatrix value, the product V/^m 
decreases rapidly with increasing j, and the resonances 
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with small values of I and m only can affect significantly a 
trajectory. Thus, if Tq > T{/i.es), only the superior term 
with I — 1 should be taken into account. Neglecting the 
higher-order resonances, we can describe the motion in 
the vicinity of /res in the pendulum approximation 5j. 
Leaving the resonant term only, we can rewrite Eas. p()|l 
in the form 



dJ 
'di 



el4(/ros)^mSin$, 



— =uj- mil + £ — ^„ cos $, 



(12) 



which corresponds with some simplification to the uni- 
versal Hamiltonian of nonlinear resonance 0, 0| 



H,, 



2^ 



[MY +eViUcos^, 



(13) 



where A/ — I — I^cs- Solutions of Eas. (|12(l describe os- 
cilations of the action / and phase <i> variables nearby 
elliptic fixed points of the respective resonances, the so- 
called phase oscillations. 

An angular location of the fixed points depends on a 
random phase ipm and, therefore, varies from one real- 
ization of the random function ^(t) to another. A tra- 
jectory of the map lO, being captured into a resonance, 
draws a chain-like structure in the phase space. In the 
regime of stable motion neighboring chains are far enough 
from each other, and the space between them is filled by 
non-resonant stable trajectories. The width of the reso- 
nance in terms of the frequency can be estimated from 
the Hamiltonian p3(l as follows: 



Alu 



1 1 ^-^max 



2J £\uj'j\Vi^r, 



(14) 



where uj'j = duj/dl is related to nonlinearity of oscilla- 
tions and determined by the form of the unperturbed 
potential. The distance between the neighboring reso- 
nances of the successive orders is 5lo{Tq) — 27r/To. If the 
criterion of Chirikov, Alo/Suj ~ 1, holds, the resonances 
overlap in the phase space, and the phase oscillations 
become unstable. 

Except for some specific situations, the Chirikov's cri- 
terion provides a sufficient condition for emergence of 
global chaos. We can estimate from the above relations 
the time that is required for overlapping all the reso- 
nances in the phase space 



(15) 



where /min is a minimal resonant action corresponding to 
the deepest resonance, and the amplitudes are sup- 
posed to be independent on time. Nevertheless, some 
stable regions in the phase space can survive at T > Tc 
in the form of islands submerged into a stochastic sea. 
The time horizon may be treated as a time of mixing 
in the phase space. With a given spectrum of noise and 



large enough values of Tq, difference between the ampli- 
tudes for different realizations of the noise becomes 
negligible. Therefore, the time of phase correlations Tc is 
independent on the realization chosen and may be treated 
as a characteristic quantity of a Hamiltonian system un- 
der consideration. 



III. CLUSTERING OF SOUND RAYS IN AN 
OCEAN WAVEGUIDE 

Low-frequency acoustic signals may propagate in the 
deep ocean to very long ranges due to existence of the un- 
derwater sound channel which acts as a waveguide con- 
fining sound waves within a restricted water volume (see 
general theory in Ref. 7| and description of field exper- 
iments in the Pacific ocean in the issue Q)- In the geo- 
metric optics approximation, the underwater sound ray 
trajectories in a two-dimensional waveguide in the deep 
ocean with the sound speed c, being a smooth function of 
depth z and range r, satisfy the canonical Hamiltonian 
equations 



dz _ dH 
dr dp ' 

with the Hamiltonian 



dp 
dr 



dH 



H = —^rfi{z, r) — 



(16) 



(17) 



where n{z,r) ~ co/c(z,r) is the refractive index, cq is 
a reference sound speed, p — nsin0 is the analog to 
mechanical momentum, and is a ray grazing angle. In 
the paraxial approximation (that takes into account rays 
launched at comparatively small grazing angles that can 
propagate in the ocean at long distance without touching 
the lossy bottom), the Hamiltonian can be written in a 
simple form as a sum of the range-independent and range- 
dependent parts 



with 



i?o = -1 + ^ 



H = Hn 



Ac(z) 

Co 



iJi(r) 



Hi{r) = 



dc{z, r) 



Co 



(18) 



(19) 



where Ac(z) — c{z) — cq describes variations of the sound 
speed along the waveguide and p ~ tan (p. The term Hi 
describes variations of the sound speed caused, mainly, 
by internal waves in ocean. 

In a simplified model of perturbation with 6c being 
a periodic function of the range r, extremal sensitivity 
of ray trajectories to initial conditions — ray dynamical 
chaos — has been fou nd [Tol [lTI| . In a number of recent 
publications 0, 0, 0, Il4j|. it has been realized that ray 
chaos should play an important role in interpretating the 
long-range field experiments 0, 0| • More realistic 
models should include not only horizontal but vertical 
modes of internal waves as well. Moreover, the internal 
waves perturbation is stochastic in the real ocean. 
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We use a background sound-speed profile introduced in 
Ref. to model sound propagation through the deep 



ocean 



c(z) = CO (^1 - y (m - e-''^)(e~"^ - 7) j , < z < h, 

(20) 

where 7 = 6 11, a, and 5 are adjusting parameters, 
h is the maximal depth, and cq — cih). We represent 
the perturbation of the background sound-speed profile 
along the waveguide in the following form: 



5c{z, r) 



Co 



eV{z)az,r), 



(21) 



where e <C 1, V{z) describes sound-speed fluctuations 
along the vertical coordinate z, and ^(z, r) is a fluctuation 
field with {S) = and (^^) — 1/2. Let us assume that the 
stochastic perturbation can be written in the factorized 
form 



az,r)^X{z)Y{r). 



(22) 



Our aim is to find coherent clusters of sound rays that 
are stable up to a distance d from a source. After making 
the canonical transformation from a ray variables (z, p) to 
the action-angle ones (/, d), we use the effective Poincare 
map © which now takes the form 



h+i=I{d- 

Let us expand Y{r) in the Fourier series 

00 

m— — 00 



(23) 



(24) 



where kr = 2'K/d, and represent X{z) as a sum of vertical 
modes 



(25) 



n— — 00 



where xiz) = e^^l^ — e^''/^, and B is a thermocline 
depth. Substituting (|24|l and (|25|l into the equations of 
motion 



dl 

dr 



Cq ' dr " ' Co dl 
one gets the resonance conditions 

-z/B 



luj — mkj. ± mrp 
Ilu + mkr ± m:p 



B 
B 



= 0, 
-0, 



(27) 



where I, m, and n take non-negative values only. 

As it is seen from (|27|l , vertical oscillations of the speed 
of sound give an additional term in the resonance con- 
ditions proportional to p and produce a special type of 



nonlinear resonances which can break regions of stability 
in the phase space 0, 0. So the question is: could 
coherent ray clusters survive under a combined horizon- 
tal and vertical stochastic perturbation? The horizontal 
perturbation Y(r) is modelled as a sum of a large (of 
the order of one thousand) number of harmonics with 
random phases and the wave numbers distributed in the 
range from /cmin = 27r/100 km~^ to fcmax = 27r km"^. 
The spectral density of the perturbation decays as k~'^. 
For more information about modelling noise see Ref. 0]- 
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FIG. 1: The map H23|l with d = 100 km revealing coher- 
ent ray clusters surviving in a model underwater waveguide 
under a noisy-like internal-wave perturbation 12211 with a 
single vertical mode 12811 . The polar action-angle variables 
~ {I / Is) cos6 and ly — {I / Is) sin6 are in units of the 
separatrix action Is- 



Let the vertical oscillations contain a single mode only 



X{z) — Xi sin TTx{z), 



(28) 



and the other parameters are as follows: e — 0.0025, 
V{z) = (2z/B)e-22/s^ and B = 1 km. We start with 
a large number of initial values of the ray action-angle 
variables in different regions of the phase space {Io,9o) 
and compute their values at the distance d = 100 km 
from a sound source with the respective equations of 
ray motion with a chosen realization of the noisy-like 
horizontal internal-wave perturbation (|21|l and the sin- 
gle vertical mode (|28f) . Then we compute the ray vari- 
ables with the same realization of noise at the same 
distance d = 100 km, find (/i,6'i), and so on. The 
result of computing of the Poincare map in the coor- 
dinates Ix = (///s)cosi? and ly = (///s)sin^ with 
Is = I{H = —1), being the separatrix value of the 
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FIG. 2: The same as in Fig. 1 but with five vertical modes 
lUnj, d = 100 km. 



action, is shown in Fig. ^ The map reveals numerous 
island-like regions of stability in the phase space corre- 
sponding to the resonance condition mD(I) ~ Id, where 
D{I) is a cycle length of a ray with a given action /. It 
should be noted that the outer chain with four islands 
corresponds to the same 1:1 resonance. Under a multi- 
frequency excitation, the system, being captured into the 
resonance mD = Id, is, simultaneously, in all the reso- 
nances jmD = jld, where j is an integer. Let us consider 
now a more complicated model with five vertical modes 
of internal waves 

5 

X{z)^J2XnSm{mTx{z) + (l)n), (29) 

n=l 

where X„ cx and (fin are random phases. The re- 

spective map H23f) is shown in Fig. [3 with d = 100 km. 
It follows from comparing Figs. and |3 that complica- 
tion in the vertical perturbation effects, mainly, the most 
steep rays corresponding to the outer resonance islands. 
The most outer island in Fig. |21 is separated from the 
other domains of stability by a layer of escaping rays 
(by escaping rays we mean the rays quiting by different 
reasons the sound waveguide during their propagation). 
With increasing d, the area of islands decreases but dis- 
tinct islands of stability remain visible in the phase space 
at a rather long distance d — 500 km (Fig. |3J) even with 
five vertical modes of internal waves to be included in the 
model. 

The practical question is: in which way the coherent 
ray clusters could be observed in real field experiments? 
What is measured in the field experiments is time ar- 
rivals of sound signals and the depths of their arrivals 
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FIG. 3: The same as in Fig. 2 but with d = 500 km. 



at a given range R with the help of a vertical receiving 
array of hydrophones 0, 0|. We compute the so- 
called timefront at i? = 1000 km solving not the map 
(|23|l but the equations of motion l|16|) with the stochas- 
tic perturbation (|24|l and (|29|l . The prominent strips 
in Fig. 2] belong, mainly, to coherent ray clusters which 
survive even at the distance of 1000 km. It should be em- 
phasized that similar strips have been found in the field 
experiments [sl flslIT^ . 




t,s 



FIG. 4: Timefront at the range 1000 km under the stochastic 
perturbation 1241 1 and 1291 1: ray depth z vs ray travel time t. 
Strips indicate coherent ray clusters. 



In conclusion of this section we show manifestations of 
coherent ray clusters in the so-called t — Po plot. Fig. \S[ 
demonstrates the ray travel time t as a function of the 
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FIG. 5; Ray travel time t vs starting momentum po at the 
range 1000 km. Shelves indicate coherent ray clusters. 



starting ray momentum po a-t the distance R = 1000 km 
from a sound source for the model with purely horizontal 
random perturbation, i. e. Ea. l22|) with X{z) — 1. The 
main feature of this dependence is "shelves" , i. e. more 
or less flat segments in the t — po plot. Each "shelf" 
corresponds to a coherent cluster. The "shelves" are dis- 
tributed chaotically over the range of the starting mo- 
menta, and their positions depend on a specific realiza- 
tion of the random perturbation. Such "shelves" have 
been found with different kinds and realizations of the 
internal-wave induced random perturbation proving that 
appearing of coherent clusters is a common feature in 
the sound-waveguide propagation through a fluctuating 
ocean. 



IV. CLUSTERING OF PASSIVE PARTICLES IN 
A TWO-DIMENSIONAL FLOW 



In this section we consider clustering in a simple two- 
dimensional flow of an ideal fluid with the dimensionless 
streamfunction 



* = In 



p,x + ex ^(r) 



(30) 



whose simplified version with ^ (r) being a harmonic func- 
tion has been introduced in Ref. j20i to model advection 
of passive particles in a flow with a flxed point vortex 
(the flrst term in ^'M)^ ). a stationary current along the 
axis y with the normalized velocity fi (the second term) , 
and a harmonically alternating current with the normal- 
ized velocity e. In physical oceanography the stream- 
function H30|l is a simple kinematic model of mixing and 
transport in the flow around a topographical eddy over 
a seamount in the ocean randomly perturbed by wind in 
the surface layer or by small-scale turbulence. It is well 
known that the Hamiltonian equations of motion for par- 
ticles (tracers) with Cartesian coordinates x and y in an 



incompressible two-dimensional flow are written as 

. _ _ y 

dy 

X 

V 



X + y 
+ li + e£,{T). 



(31) 



dx x"^ + 

The configuration space of advection particles is the 
5 
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FIG. 6: The map @ with To = Gtt revealing clusters of ad- 
vected particles in a random velocity field with the amplitude 
e = 0.01 and the frequency band uj £ [2:5]. The unperturbed 
frequency of rotation of particles around the point vortex ujq 
and the polar angle tp/TT are in arbitrary units. 



phase space of dynamical system (|30ll . and it is possi- 
ble to see noisy-induced clusters of particles in the flow 
(|30|l by a naked eye. We model the noise as a random 
function ^(t) consisting of a large number of harmonics 



1 ^ 

Fij) = ; siii(w„r + (fin) 

with random phases (p„ and the frequencies 



N 



0,1, 



(32) 



(33) 



distributed equally in the range [ujb '■ i-^e] between the 
lowest (uJb) and highest (wg) frequencies. It is possible, 
in principle, to model the spectrum of an arbitrary form 
with the help of the series H32|l introducing an amplitude 
factor depending on the frequency. Accordingly to the 
central limit theorem, ^(t) has a Gaussian distribution 
with zero expectation value because it is a large number 
of independent random variables. The factor + 1 

provides the variance to be equal to 1/2. 

Advection of passive particles under a single- mode per- 
turbation without any noise has been considered it detail 
in Rcfs.pll.l23| as a scattering problem: particles enter a 
mixing zone around the vortex with the incoming current 
and escape from it with an outcoming flow. The scatter- 
ing has been shown to be chaotic in the sense that there 
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FIG. 7: The averaged (over To) energy of particles (E) vs their 
unperturbed frequency luq: (a) To = 6tt and (b) To = SOvr. 
Shelves indicate clusters of particles. Parameters of noise are 
given in caption to Fig. |S| 



is a nonattracting invariant set, consisting of unstable 
periodic, aperiodic, and marginally unstable orbits, that 
determines scattering and trapping of particles from the 
incoming flow. The scattering has been shown to be frac- 
tal in the sense that scattering functions (namely, the 
dependence of the trapping time for particles and of the 
number of times they wind around the vortex on the ini- 
tial particle's coordinates) are singular on a Cantor-like 
set of initial positions. 

Now we are interested in what happens if that deter- 
ministic flow is perturbed by a weak noise. It should be 
noted, first of all, that fractality in the scattering func- 
tions does not disappear under noisy excitation, but the 
respective fractals have been shown to have more com- 
plicated hierarchical structure [l^ . In order to find 
noisy-induced clusters of particles we use the effective 
Poincare map of the type of © with different values of 
the time interval Tq, the amplitude of noise e, and the 
frequency range of noise [cjf, : LUe]. 

In Fig. we demonstrate the clusters of the map of 
the type © with Tq = 6n in the plane of the frequency 



luq of rotation of particles around the vortex in the un- 
perturbed system with e — and the polar angle with 
tamp — y/x. A large number of particles advected by 
the flow H3U|) with = 0.5 and a chosen realization of the 
weak noise with the amplitude e ~ 0.01 and the medium- 
frequency range w G [2 : 5] has been chosen for simula- 
tion. 

Indirect manifestations of clusters are seen in Fig. 
as shelf-like segments in the dependence of the averaged 
(over To) energy of tracers {E) on their unperturbed fre- 
quency Wo at initial positions. With increasing the in- 
terval of mapping Tq, noisy-induced clusters reduce, of 
course, in their size, but even with Tq = SOtt, which is 
larger than the average period of perturbation almost by 
two orders of magnitude, the prominent shelves corre- 
sponding to small-size clusters are still seen in Fig. 03- 
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FIG. 8: The map Q with To = 2.37r revealing clusters of 
particles in a random velocity field with the amplitude e = 0.1 
and the frequency band a; G [9 : 10]. (a) general view in the 
phase plane {x, y), (b) and (c) magnifications of the respective 
clusters in (a) . The dotted line is an unperturbed separatrix. 



8 




cOq (a.u.) 

FIG. 9: The averaged (over To) energy {E) of particles vs their 
unperturbed frequency ojq: (a) To = 2.37r and (b) To = IOtt. 
Shelves indicate clusters of particles. Parameters of noise are 
given in caption to Fig. |H| 



In order to test the effectiveness of the map of the 
type © in "catching" the noisy-induced clusters, we have 
carried out a series of numerical experiments with /i = 0.5 
and the noise of the moderate amplitude e = 0.1 and the 
high-frequency range a; G [9 : 10]. Fig. demonstrates 
the result of mapping with Tq — 2.37r on the configuration 
plane {x,y), and Figs. and c are magnifications of 
the respective clusters in Fig. Prominent chains of 
the islands of stability can survive under a rather strong 
noisy perturbation for a comparatively long time. Fig. |51 
demonstrates shelf-like manifestations of noisy-induced 
clusters with Tq — 2.3tt and Tq = IOtt. 

To give a more direct manifestation of coherent clus- 
tering in a random flow that could be observed in real 
laboratory experiments with dye, we compare the evolu- 
tion of patches of particles chosen in the regions of sta- 
bility and instability in the configuration space. In the 
upper panel of Fig. ^] we show the evolution of a coher- 
ent cluster corresponding to the small black point (the 
region of stability for a given realization of noise) with 
coordinates x ~ —0.55 and y ~ 0.1 in Fig. ISji. More or 
less compact evolution of this patch with lO"' particles 
goes on up to, at least, t = 6tt. For comparison, the 
lower panel of Fig. llOl demonstrates on the same time in- 
terval the evolution of the patch with the same number of 



particles chosen initially close to the coherent cluster at 
T — 0. The respective initial patch is deformed strongly 
at T = Gtt. 

Clustering in the ocean is a common feature that can 
be seen, for example, in satellite images of the ocean 
surface. In resent years new observational tools — quasi- 
Lagrangian current following fioats and drifters — have 
been used to observe velocity field in the ocean at differ- 
ent levels of depth (for a review see H^l). In connection 
with our numerical observation of coherent clusters of 
passive particles in the simple kinematic ocean model, 
we would like to pay attention to the results of the SO- 
FAR fioats program in the POLYMODE experiment in 
the North Atlantic 251- Up to forty neutrally buoyant 
floats at 700 m and 1300 m were used to provide a quasi- 
Lagrangian description of the structure and evolution of 
the mesoscale eddy field. Those floats can be considered 
as quasi-passive tracers in a weak noise environment. A 
large number of the deep floats revealed remarkably co- 
herent motion over a two-month period. 



V. CONCLUSION 

We proposed and justified a rather general method to 
find regions of stability in the phase space of oscillatory- 
like Hamiltonian systems driven by a weak noise with 
an arbitrary spectrum. Physical manifestations of these 
regions of stability, the so-called coherent clusters, have 
been demonstrated with two models in ocean physics. 
We have found coherent ray clusters in the model of 
underwater-waveguide sound proparagition through a 
randomly-fiuctuating ocean and coherent clusters of pas- 
sive particles in the two-dimensional flow modelling kine- 
matically advection around a topographic eddy in the 
ocean. 

In conclusion we would like to add two remarks about 
the properties of the effective Poincare map ©■ Firstly, 
the duration of the temporal interval Tq in construct- 
ing the map can be chosen arbitrarily. Thus, if ^{t) is 
a stationary random process, the regions of stability in 
the phase space exist at any time moment. Secondly, 
the map © enables to prove definitely the existence of 
some regions of stability but not all of them. Really, the 
map, constructed with any given value of Tq, can reveal 
only those stable sets which correspond to the phase os- 
cillations nearby the fixed points of the map. However, 
there exist another regions of stability looking as chaotic 
ones on the map. The topology of the map changes with 
varying the mapping time Tq, and some regions in the 
phase space, which look as pseudochaotic on the map at 
To = Ti , become stable at Tq = T2 > Ti . The total area 
of the regions of stability, survived under a weak noise 
at i = To, can be estimated as an area of superposition 
of all the stable sets detected by the map with the 
mapping step varying from Tq to Tc- 
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FIG. 10: The compact evolution ol a coherent cluster of particles (the upper panel) and deformation of an unstable patch of 
particles (the lower panel) in a random velocity field. Parameters of noise are given in caption to Fig. |H| 
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